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Abstract 

In this letter we construct the Stiickelberg holographic superconductor with Weyl corrections. 
Under such corrections, the Weyl coupling parameter 7 plays an important role in the order of 
phase transitions and the critical exponents of second order phase transitions. So do the model 
parameters c a , a and C4. Moreover, we show that the Weyl coupling parameter 7 and the model 
parameters c a , a, C4 which together control the size and strength of the conductivity coherence 
peak and the ratio of gap frequency over critical temperature co g /T c . 
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I. INTRODUCTION 



As is well known, there exist lots of questions in strongly coupled field theories. Facing 
these problems, many physicists resorted to the AdS/CFT correspondence (lHj|, a equiva- 
lence between a d- dimensional strongly coupled gauge theory on the boundary and a (d+ 1)- 
dimensional weakly coupled dual gravitational description in the bulk, to deal with some 
unsolved problems. Nowadays, as an application of the correspondence, the AdS/Condensed 
Matter duality (AdS/CMT) is also being developed. One remarkable achievement of AdS/ 
CMT is the establishment of the holographic superconductors 5j. It originated from the 
discovery of Abelian Higgs mechanism, the U(l) gauge symmetry is spontaneously broken 
due to the existence of a black hole p, ]7|. For a general review of the subject, please see 

St 

Of course, this model also displays some interesting properties. For example, there is a 
second order phase transition which is the standard behavior shown by Landau-Ginzburg 
theory. Another important feature is that the ratio of gap frequency over critical temperature 
oj g /T c Rj 8, which is much larger than the value of 3.5 given by the weakly coupled BCS 
theory, but it is quite similar to that in the high-T c superconductors experiment 13] . 

Recently, Franco et al. has proposed a phenomenological model, in which the spontaneous 
breaking of a global U(l) symmetry occurs via the Stiickelberg mechanism 14|-|l6l|. They 
regard this model as the Stiickelberg holographic superconductor which has lots of features. 
One main characteristic of this phenomenological model is that it can provide a description 
of a large group of phase transitions, including the first order phase transition and the second 
order phase transition with both mean and non-mean field behavior. Many extending works 
on the Stiickelberg holographic superconductor have been explored, such as the Stiickelberg 
holographic superconductor with Gauss-Bonnet corrections, the Stiickelberg holographic su- 

and the Stiickelberg holographic superconductors 



perconductor with backreactions 
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in constant external magnetic field 19j. In Ref. [l7j, they find that different values of GB 
parameter and model parameters can determine both the order of phase transitions and the 
critical exponents of second order phase transitions. What's more, the size and strength of 
the conductivity coherence peak can be controlled. 

In this paper, we pay attention to the effects of the Weyl corrections, another higher 
derivative corrections on the Stiickelberg holographic superconductor. The Weyl corrections 
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involve a coupling between the Maxwell field and the bulk Weyl tensor. In Ref. [201 ] . 
the authors calculated the conductivity and charge diffusion with Weyl corrections, and 
showed that corrections broke the universal relation with the U(l) central charge observed 
at leading order. In addition, the holographic superconductors with Weyl corrections are 



also explored in Ref.[2l|. They found that for large Weyl parameter 7, the ratio of gap 
frequency over critical temperature oj 9 /T c is far less than the value 8 which is found in 
the standard version holographic superconductors 0, l^j, while for small 7, the value of the 
ratio is larger than 8. Therefore, it is worth exploring how the picture of the Stuckelberg 
holographic superconductor is modified when considering Weyl corrections. 

Our paper is organized as follows. In section II, we construct Stuckelberg holographic 
superconductor with Weyl corrections and present numerical results from the Stuckelberg 
holographic superconductor. We specially discuss the order of phase transitions and the 
critical exponents of second order phase transitions. Then we investigate the conductivity 
numerically in section III. Conclusions and discussions are followed in section IV. 

II. THE STUCKELBERG HOLOGRAPHIC SUPERCONDUCTOR WITH WEYL 
CORRECTIONS 

A. The Stuckelberg holographic superconductor with Weyl corrections 

The general Stuckelberg holographic superconductor contains the gravitational field, the 
U(l) gauge field and the scalar field, the latter two of which are coupled via a general 
Stuckelberg mechanics Il5 |. The action in this system is 



S= l^zV^lR+^j^-C), (1) 



with 



where J 7 is a general function of ^> 

J-(v&) = ^ 2 + c Jf a + c 4 ^ 4 , (3) 

in which ca. c a , a are real numbers and 3 < a < 4. When c a = c 4 = 0, it reduces to the 
model of [5|. In addition, this theory has the local gauge symmetry — > + c^f^x^) 



and p — > p + where we can use the gauge freedom to fix p = 0. Furthermore, we take 

the ansatz $ = W. Therefore, the Stiickelberg Lagrangian ([2]) can be re-expressed as 

L = + JL ^~ + -y* 2 + 2-?W ( 4 ) 

The function J 7 becomes 

JT(^) = ^ 2 + Ca ^« + c 4 # 4 . (5) 

In this paper, we will construct a 3 + 1 dimension Stiickelberg holographic superconductor 



21] . In this case, the matter Lagrangian density L is replaced 



with Weyl corrections as Ref. 
by 

C w = - {F^F, U - AnC^F^) + + — ^> 2 + -.F(¥) V*", (6) 

where 7 is a (dimensionful) constant and C M „ pCT is the Weyl tensor. The range of the param- 



eter 7 [20( is 



L 2 L2 

< 7 < — , (7) 

16 ' 24' v ; 

where the upper bound results from the existence of an additional singular point at 7 = |j 

and the lower bound is due to the causality constraints. 

We exclusively work in the probe approximation in this letter, where the gravity sector 

is effectively decoupled from the matter sector. So we take the background to be a planar 

AdS-Schwarzschild black hole: 

ds 2 = ^-(-f(r)dt 2 + dx t dx l ) + Sr^dr 2 , (8) 



where 

4 

J 

with Tn > for the horizon and r — y 00 for boundary of the bulk, the Hawking temperature 
of this black hole is 

r = £, (io) 

which corresponds to the temperature of the dual field theory. For convenience, in the 
following analysis, we always set L = 1. Let's begin with the following Ansatz: 

tt = tt(r), At = $(r), (11) 



then apply it to the Euler-Lagrange equation, we derive the equations of motion from the 
Weyl corrections matter Lagrangian density ([6]) 

-^)*"+(^)*'-^0, (!3) 

where the prime represents the derivative with respect to r and the dot in T represents the 
derivative with respect to In addition, we have used the nonzero components of Weyl 
tensor C^p^ in AdS^, 

Coioj = f(r)r H 5ij, Coo = — Cirjr — ~ Cijki — i" H 5ik5ji. (14) 

For concreteness, we fix the mass of the scalar field to m 2 = —3, which is above the 
Breitenlohner-Freedman bound 



23] and does not induce an instability With this mass, 
the asymptotic behavior of the scalar fields ^> and the scalar potential $ close to the AdS 
boundary (r — > oo) take the form 

• = £ + £ (") 

* = f-£. (16) 

where /i and p are the chemical potential and charge density in the dual field theory, respec- 
tively. For the value of m 2 we have mentioned before, only the ^3 term is normalized. Here 
we set = 0, so the expectation value of the scalar operator O3 in the dual field theory is 
related to ^3 by (O3) = ^3. Furthermore, at the horizon r = rjj, the regularity implies 

nrn) = ~*(r H ), (17) 

Hr H )=0. (18) 

To investigate the behavior of phase transitions of this system, we solve equations of 
motion (fl2l) and (fi3|) combined with the boundary conditions (jl"5l) . (fl"6l) . (TTTT) and (TIB"]) 
numerically by using a shooting method. 
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B. Phase transitions 



One main characteristic of the Stuckelberg holographic superconductor is that it enables 
us to describe various types of phase transitions, from first to second order. As for the 
second one, it is also able to illustrate those with the mean and non-mean field behaviors. 
When involving the Weyl corrections, how the phase transition depends on the coefficients 
c a , a, c 4 and Weyl coupling parameter 7 is deeply investigated. In order to see the role each 
parameter plays, we observe the following three cases respectively. 



1. The case JF{^) = V 2 + c 4 ^ 4 

Firstly, we focus on the special case J r (^ r ) = \1/ 2 + c 4 \I/ 4 to see how c 4 and 7 control the 
styles of the phase transition. Solving the equations of motion ( Tl2|) and ([131) numerically, we 
have the results itemized as follow: over the legal range of 7, the phase transition is second 
order when < c 4 < 0.7 and is first order for c 4 > 1.3 (Fig{T]). However, if 0.8 < c 4 < 1.2, 
the phase transition varies from second order to first order, and first order ones are prone 
to appear in smaller 7, while the second order ones usually emerge in larger 7 (FigJT]). 
Furthermore, by fixing the data, we find that for the second order phase transition, while 
for T ~ T c , there is a relationship between C 3 and T c , that is, O3 ~ (T c — T)' 3 , = 1/2, 
which has nothing to do with 7. 



2. The case F{^>) = 4< 2 + c Q ^ a 

Now, in order to investigate the influence of c a , a and 7 on the phase transition, we set 
c 4 = 0. In Figf2[ we present the condensate around the critical region for chosen values of 7, 
c a and different values of a. Fig f5] informs us that for chosen values of 7 and c a , first order 
phase transitions often gather in the region of smaller a, whereas second order ones usually 
lie in the region of larger a. Moreover, for fixed 7 and a, the styles of phase transition exhibit 
some kinds of similar distribution, but what is different is that first order ones appear with 
larger c a , while second order ones happen when c a is small. To sum up, more often than 
not, the first order phase transition locates in the region of larger c a or smaller a, which is 



consistent with the result in Gauss-Bonnet gravity 17| • Then, we will specially concentrate 



on the effect of the Weyl parameter 7 on the phase transition. From Figj2l we find that 
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FIG. 1: The condensate around the critical region for chosen values of C4 and different values of 
7. To is the temperature when condensate and second order phase transition appear, To equals to 
the critical temperature T c . The orange, black, magenta, green, red and blue curves correspond to 
7 = -0.06, -0.04, -0.02, 0, 0.02, 0.04, respectively. 

for chosen values of a and c a , first order phase transitions are distributed in small 7, while 
second order ones occupy the larger 7 region. In light of discussions above, one can conclude 
that the first order phase transition of the system is liable to emerge when Weyl parameter 
7 is small. 

3. The case JF{^) = + 

In the above two cases, we have shown that different model parameters can lead to either 
first or second order phase transitions. It is natural to ask whether the change of model 
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FIG. 2: The condensate around the critical region for chosen values of 7, c a and different values 
of a. The red, blue and green curves correspond to a = 3, 3.5, 4, respectively. 



parameters would exert some influence on the critical exponent, rende ring it different from 
the prediction of the mean filed. In reality, the answer is yes. Ref. 15] has shown that 
for c a < 0, the critical exponent is not the same as the result the mean field holds when 
model parameters are varied. When involving the Weyl corrections, we also focus on the 
special case c a = — 1 and C4 = 1/2 as in Ref. 15J. From Fig|3j we find that the critical 
exponent /3 is independent of the Weyl parameter 7 but only dependent on a. Adopting 
the non-linear fitting method, we can get the critical exponent /3 for different a. Therefore, 
the relation /3 = (a — 2)~ 1 found in Ref. 15[ is still preserved. These results show that the 



critical exponent /3 depends only on the model parameter a but is independent of the Weyl 
parameter 7. Actually, according to other parameters, such as the scalar mass, the 

dimensions of space-times and the background geometry, have no influence on the value of 
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FIG. 3: The critical exponent (3 (the slope of the curves) is independent on the the Weyl parameter 
7 and only dependent on a. The vertical axis is the log of the Condensate (O3) for the case 
near the critical temperature T c and the horizontal axis is the log of 1 — T/T c . 
The blue, red and green curves correspond to 7 = —0.06, 0.02, 0.04, respectively. 

the critical exponent (5. Thus, we can say that in this certain case Ffif) = + |^ 4 , the 
critical exponent is decided by parameter a. 



a 


3 


3.25 


3.5 


4 


P 


1 


0.79 


0.66 


0.5 



TABLE I: The critical exponent /3 for the case J 7 ^) = —^> a + with different a. 

Next, we will focus on the dependence of the condensate on a and 7. First, we present 
the condensate as a function of temperature in FigJH From these figures above, it is easy 
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FIG. 4: The condensate around the critical region for chosen values of 7, c a and different values 
of a. The red, blue and green curves correspond to a = 3, 3.5, 4, respectively. 

to see that the given different values of a and 7, the critical temperature T c are unrelated 
to a and only rely on 7. In the meantime, the critical temperature T c for different Weyl 
parameter 7 is summarized in Table [Til Observing this table, we claim that the critical 
temperatures T c is higher and higher as the Weyl parameter 7 changes from —0.06 to 0.04, 



which is in agreement with Ref. [2l|. In another word, when Weyl parameter 7 runs toward 
the negative direction, it is much harder for the scalar hair to form. However, if 7 turns 
around, the formation of the scalar hair would be facilitated. In fact, for the second order 
phase transition in all cases above, the critical temperature T c only depends on the Weyl 
parameter and the form of J 7 ^) is not the determinant of T c . 
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7 


-0.06 


-0.02 


0.02 


0.04 


T c 


0.170P 1 / 3 


0.185P 1 / 3 


0.219P 1 / 3 


0.30V/ 3 



TABLE II: The critical temperature T c for the case J~(^) = —^/ a + with different 7. 



III. ELECTRICAL CONDUCTIVITY 



In this section, we study the electrical conductivity a of the dual field theory. We 
switch on a vector perturbation of radial symmetry and time dependence in the x-direction, 
SA x (t, r) = A x (r)e~ tu3t dx 1 , and study the linear response of the system. With the help of 
the Euler-Lagrange equation, the equation of motion for the perturbation of vector potential 
SA x (t,r) is obtained in the following 

(H)+^(f 1 

87?^ 



:i + 



A" 



+ 1 + 



—A 



r 4 f 



r 2f x 



r 

-- 0. 



AL 



(19) 



At the horizon, the ingoing boundary conditions due to the requirement of causal behavior 



24J must be imposed. In our case, this is 

Mr) ~ (r 2 f(r))~ 1 ^, 
Around the AdS boundary r — > 00, A x (r) is expanded as 

Mr 



(20) 



A(0) | A^ | A^WlogAr | 



(21) 



r 2 r z 

where A@\ A^ and A are integration constants. The logarithmic term suffers from an 
ambiguity scale A, which leads to a divergence in the Green's function. But as explained in 



22j, this can be removed by introducing an appropriate counter-term in the action. From 
linear response theory, the conductivity a is related to the retarded current Green's function 
G R as follow, 



a (u) = ^G R (u,k = 0), 



(22) 



1 for simplicity, we will work at zero momentum (k = 0). 
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FIG. 5: Conductivity for superconductors with different values of a, c a and 7 at the same temper- 
ature T/T c ~ 0.4. The solid line represents the real part, and dotted line the imaginary part. 
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FIG. 6: Re(a) as a function of uj in units of the condensate with different values of a, c a and 7 at 
the same temperature T/T c « 0.7. 

where the retarded current Green's function G R can be evaluated by using the standard 
AdS/CFT technique as 2 

G» = - hm r*f{r) *ff>"\ (23) 

Substituting the solution (12 ip into the above equation, one can easily obtain the explicit 
expression of the retarded Green's function 



A (2) 



0J 



2 



Therefore, the conductivity is given by 

° M = -^W)+r (25) 

In order to read off and A^ 2 \ we must solve Eq.( jl"9l) . Eqs.( fT2l and (fT3l) numerically 
with the boundary condition ( l20l) . For definiteness, we set C4 = and explore the effect of a, 
c a and 7 on the conductivity. The numerical results of the frequency dependent conductivity 
are illustrated in FIGj5]and[6j The main characteristics are summarized as follows: 

1. The imaginary part has a pole at u — 0. According to the Kramers-Kronig relation, 
this indicates that the real part contains a delta function. In addition, from FIGJ3J, 
it is easy to find that there exists a gap in the conductivity, which rises quickly near 



2 In Rcf.[2l[, the authors have given a detailed derivation of the retarded Green's function. A similar 
calculation can be also found in Ref-f^. They find that Wcyl term has no explicit effect on the retarded 
Green's function and its expression is the same as that in Einstein theory for a standard Maxwell field. 
Therefore, in the following calculation, we can safely use the Eq. 
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the gap frequency . g . They aretwo main common features of the standard version 
holographic superconductors |5l. I22|. 

2. The ratio of gap frequency over critical temperature ou g /T c is unstable and running 
with the model parameter a, c a and 7 (FIGJ5]). For fixed values of a and c a , the ratio 
cu g /T c increases with the fall of the Weyl parameter 7. Especially, for 7 = 0.04, the 
value of ratio QjJTr. y 6, far less than 8 found in the standard version holographic 



superconductors 



22] . which is sharply different from the previous results. Therefore, 
compared with the former models, the value of our system is closer to the weakly 
coupled BCS value of 3.5. Moreover, if we fix values of Weyl parameter 7, the ration 
a; g /T c will increase with the augmentation of model parameter a or c a . 

3. In Ref. 



2l|, they have found that an extra spike appears inside the gap when 7 < 0. 
But for this model, the new things are not only when 7 < but also when a or c a are 
large, the spike will emerge, too. Moreover, FIGJ5] implies that under the joint effects 
of 7, a and c a , that is, when 7 < and a and c a are large, there will be more spikes 
and each of them appears to be sharper and narrower than those under the effect of 
one single model parameter. 

4. The parameter 7, a and c a together control the strength of fluctuations in this system. 
For fixed value of 7, the coherence peak gradually becomes stronger and narrower 
with the increase of a or c„ (FIGJ61), which is in agreement with the result found in 

nn 

Refs.|15l. 117] . For fixed value of a and c a , the coherence peak becomes stronger and 
narrower with the fall of the Weyl parameter 7 and the fluctuations will be suppressed 
completely at larger value of 7 = 0.04 (FIG]6]). 



IV. CONCLUSIONS AND DISCUSSION 



In this paper, we have constructed the Stiickelberg holographic superconductor with the 
Weyl corrections and obtained the numerical solutions of this model. We examine carefully 
the effects of model parameters, c a , a, C4 and the Weyl coupling parameter 7 on the phase 
transitions. For the second order phase transition, the critical temperature T c only depends 
on the Weyl parameter and the form of J 7 ^) is not the determinant of T c . For fixed model 
parameter c a , a and C4, the critical temperature will be increasingly high as we amplify the 
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parameter 7. Therefore, the condensation gets harder when 7 < 0. But for 7 > 0, the 



2l|. 



tendency is just opposite. It is in line with the results found in Ref. 

We specially focus on how the Weyl coupling parameter 7, together with the model 
parameters c a , a, C4, exerts an effect on the order of phase transitions and the critical 
exponents of second order phase transitions. We find that the Weyl coupling parameter 
7 also changes the order of the phase transition besides the model parameters c a , a, C4. 
Furthermore, we conclude that for some definite form of Ffif), the critical exponents /3 is 
only dependent on the model parameter but independent of the Weyl coupling parameter 7 
in second order phase transitions. 

Finally, we calculate the conductivity of the Stiickelberg holographic superconductor 
with the Weyl corrections numerically and obtain that for fixed model parameters c a , a and 
C4, the ratio u) g /T c is unstable and becomes smaller when the Weyl coupling parameter 7 
increases, and even falls to approximately 6, which is far less than the standard value 8. It is 



in agreement with the discovery in Ref. [21] . In addition, we also examine how the parameter 
7, a and c a control the strength of fluctuations in this model together. We observed that 
the strength of fluctuations was adjusted not only by the parameter a and c a as found in 



Refs. 



but also by the Weyl coupling parameter 7. 
On the whole, the higher curvature interactions on the gravity side broaden the class of 
the dual field theories one can holographically study. Various holographic superconductors 



model taken into account the higher curvature interactions are working intensively 26 



331 ] . It may be interesting to further explore the joint effects of various higher curvature 
interactions, including Gauss-Bonnet, quasi-topological and Weyl terms, on the holographic 
superconductors. 
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